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Abstract 

A morph between two Riemannian n-manifolds is an isotopy between 
them together with the set of all intermediate manifolds equipped with 
Riemannian metrics. We propose measures of the distortion produced by 
some classes of morphs and diffeomorphisms between two isotopic Rie- 
mannian n-manifolds and, with respect to these classes, prove the exis- 
tence of minimal distortion morphs and diffeomorphisms. In particular, 
we consider the class of time-dependent vector fields (on an open subset 
of E n+1 in which the manifolds are embedded) that generate morphs 
between two manifolds M and N via an evolution equation, define the 
bending and the morphing distortion energies for these morphs, and prove 
the existence of minimizers of the corresponding functionals in the set of 
time-dependent vector fields that generate morphs between M and TV and 
are L 2 functions from [0, 1] to the Sobolev space w£ ' 2 (fi, R n+1 ). 

1 Introduction 

1.1 Summary of Results 

Let M and TV be compact and orientable smooth Riemannian n-manifolds iso- 
metrically embedded into M n+1 . A morph between M and TV is an isotopy 
between them together with the set of all intermediate manifolds equipped with 
the Riemannian metrics inherited from Every morph or diffcomorphism 

between isotopic manifolds produces distortion via stretching and bending. We 
define functionals that measure distortion and prove the existence of minimal 
distortion morphs and diffeomorphisms. 

Let f2 C R ra+1 be an open set containing the manifolds M and TV. We define 
functionals E and £ that measure the distortion of diffeomorphisms and morphs 
respectively generated by time-dependent vector fields v : ft x [0, 1] — > R n+1 via 
the evolution equation dq/dt — v(q, t) and prove the existence of minimizers of 
E and £ in an admissible set A k P of time-dependent vector fields, which is a 



subset of the closed ball of radius P in the Hilbert space TL k of all L 2 functions 
from [0, 1] to the Sobolev space Wq' 2 (£1; R" +1 ), where k £ N. We also analyze in 
detail a concrete example of a minimal morph for the case of circles embedded 
in the plane. 

We compute the Euler-Lagrange equations for the deformation energy func- 
tional defined on Diff (M, N) and show that the radial map between a manifold 
and its rescaled version is a critical point of the deformation energy functional. 
We prove the existence of minimal deformation holomorphic diffeomorphisms of 
Ricmann surfaces. 

1.2 Background and Motivation 

A fundamental problem in Riemannian geometry and related areas is to deter- 
mine whether two diffeomorphic compact Riemannian manifolds (M, c/m) and 
(N,gpf) are isometric; that is, if there exists a diffeomorphism h : M — > N 
such that h* gN — gu — 0, where h* g^ denotes the pull-back of <?at by h. If 
no such diffeomorphism exists, it is important to know whether there exists a 
diffeomorphism that most closely resembles an isometry. This is accomplished 
by minimization of the deformation energy functional 



over the space Diff (M,N) of diffeomorphisms between M and N, where u>m is 
the volume form on M and || ■ || is the fiber norm on the bundle of all (0, 2)-tensor 
fields on M generated by the fiber metric g* M ® g* M . 

The minimization problem takes on added significance once the physical in- 
terpretation of the tensor h* g^ — g^ is recognized: it is exactly the (nonlinear) 
strain tensor corresponding to the deformation h in case gu and g^ are Rie- 
mannian metrics inherited from Euclidean space. Thus, this functional and its 
variants must occur in physical problems. Indeed, the minimal distortion prob- 
lem arises, for example, in manufacturing, computer graphics, movie making, 
and medical imaging. 

The problem of bending a sheet of metal to a desired shape using minimal 
energy has been studied (see |22j . [?]); but, algorithms for numerical approxi- 
mations are proposed and used without proving the existence of minimizers. 

An animation might require an aesthetically pleasing transformation that 
takes one image to another through intermediate shapes. Such a transformation 
is called a morph, or a metamorphosis (see |21j for a survey on morphing). 
A desirable morph might be defined as a minimizcr of a cost functional that 
measures the distortion energy. In [TS], the distortion energy of an elastically 
dcformable surface r : U x [0, T] — > K 3 , where U C K 3 is an open set, produced 
by a deformation h — r(-, T) is defined as the integral of weighted norms of the 
local coordinate representations of the strain tensor h*gx — go an d the tensor 
h* IIt — Hoi where gt and lit are the first and the second fundamental forms of 
the surface r(U,t) at the time t £ [0,T]. 
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Image matching and image registration is an important subject in medical 
imaging. Image matching is used for determining the existence of abnormali- 
ties (distortions due to underlying medical conditions) in two images, taken at 
different times, of the same organism. The problem of registration of a popula- 
tion of images to one template for the purpose of statistical analysis is another 
instance of image registration. One approach to the image matching problem is 
by minimization of a distortion functional (see [51 151 HTJl ITU] ) . 

In the following sections, we will discuss the underlying mathematical prob- 
lem of the existence of minimal distortion diffeomorphisms and morphs between 
embedded manifolds of codimension one. 

1.3 Mathematical Preliminaries, Definitions, and Results 

Given a smooth oriented n-manifold S (perhaps with boundary) isometrically 
embedded into K" +1 , we let gg, ws, and lis denote the Riemannian metric, 
volume form, and second fundamental form on S associated with this embed- 
ding. Also, we let Int S (respectively, dS) denote the interior (respectively, the 
boundary) of the manifold S. 

Definition 1.1. Let M and N be isotopic compact connected smooth n-manifolds 
(perhaps with boundary) embedded in R n+1 such that M is oriented. A C°° 
isotopy F : M x [0, 1] — > K™ +1 together with all the intermediate manifolds 
M* := F(M,t), equipped with the orientations induced by the maps /* = 
F(-,t) : M — > M l and the Riemannian metrics g t inherited from R n+1 , is called 
a (smooth) morph from M to N. 

We denote the set of all smooth (respectively, C r ) diffeomorphisms between 
manifolds M and N by Diff (M, N) (respectively, Diff (M, N)). Similarly, we 
denote the set of all smooth morphs between M and N by A4(M,N). If F is 
an isotopy, then each map F(-,t) : M — > M* induces smooth diffeomorphisms 
IntF(-,t) : IntM -> IntM* and dF(-,t) : dM -> 9M* by restriction. 

In addition, we consider morphs between manifolds M and N with different 
regularity properties. For example, we let A4 r ' ac (M, N) denote the set of all 
continuous isotopies F : M X [0, 1] — > M. n+1 between M and such that for 
each p £ M the map 1 1— > F(p, t) is absolutely continuous on [0, 1] and for each 
t £ [0,1] the function p i— > F(p,t) is a C' diffcomorphism from M onto its 
image. As in the case of smooth morphs, the diffeomorphism F(-,t) : M — > M* 
induces an orientation on the intermediate manifold M*. 

There are several choices for cost functionals that measure the distortion of 
a diffcomorphism h £ Diff (M, AT) or a morph J 1 £ A4(M, N). 

A complete theory of the existence of minimizers of cost functionals that 
measure distortion of diffeomorphisms and morphs due to change of volume is 
presented in [5]. In this case, the value of the distortion energy functional at 
a diffcomorphism h : M — * N is defined to be the square of the infinitesimal 
relative change of volume \ J(h) \ — 1 produced by h integrated over the manifold 
M, where J(h) is the Jacobian determinant of h. This functional does not 
take into account the distortion of shape produced by h, which is captured by 
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functionals and where the fiber metric || • || on the bundle of all tensor 
fields of type (0, 2) is induced by the fiber inner product g* M ® g* M (see [T5]). 

The general problem of the existence of minimizers of $ is open. The special 
case where M and N are one-dimensional is studied in [5] where, among other 
results, the functional $ is shown to have no minimizer in case M and N are 
circles with the radius of N smaller than the radius of M . Thus, a solution of 
the general problem must take into account at least some global properties of 
the metric structures of the manifolds M and N. On the other hand, we will 
prove the existence of minimizers in case M and N are Riemann spheres or 
compact Riemann surfaces of genus greater than one and the admissible set is 
HD(M, N) = {h e Diff(M,iV) : h is a holomorphic map}. More precisely, the 
following theorem will be proved in section [4] 

Theorem 1.2. (i) Let fiR : ]R 3 — > M 3 be the radial map given by hn(p) = Rp for 
some number R > 0. If M = S 2 is the 2- dimensional unit sphere isometrically 
embedded into R 3 and N = hn(M), then h G HD(M, N) is a global minimizer 
of the functional restricted to the admissible set HD(M, N), if and only if 
h = f o IirIm) where f is an isometry of N . 

(ii) Let M and N be compact Riemann surfaces. //HD(M, N) is not empty and 
the genus of M is at least two, then there exists a minimizer of the functional 
$ in HD(M,7V). 

We note that the diffeomorphism h — f o hji \ m in the latter theorem satis- 
fies the Euler-Lagrange equations for the functional <f> with its natural domain 
Diff(M, TV) (see appendix [A| . The Euler-Lagrange equations for the deforma- 
tion energy functional <!> : Diff (M, N) — » M. + are highly nonlinear and rather 



complicated (see proposition A.3l, which discouraged us from using them to 
show the existence of minimizers of the functional <I> on Diff (M, N). 

If we wish to match, in addition to the Ricmannian metrics, the embeddings 
of the manifolds M and N (to avoid, for example, zero distortion energy maps 
between a square and a round cylinder in M 3 ), we arrive at the problem of 
minimization of the functional 

2 , . i / II U* TT_. TT_ . II 2 , 



A(h):= \\h*g N -g M fu M + \\h* ll N - ll M \\ 2 u M (2) 
Jm jm 

over the space of diffeomorphisms between M and N, where IIm and IIjv are 
the second fundamental forms on the manifolds M and N. 

One of the difficulties encountered in attempts to minimize $ over Diff (M, N) 
is the lack of a complete understanding of the structure of this infinite-dimensional 
space. The natural new approach is to linearize; that is, replace Diff (M, N) with 
a subset of a linear function space. Using this approach, which already appears 
in the literature on image deformation (see [BJ |5J EHJ EH]), we define our distor- 
tion energy functionals on time-dependent vector fields that generate morphs 
(see Fig. [I). 

Let us denote the Euclidean norm of an element A € E m by \A\ or by 
|A|n.m. Let C E™ +1 be an open ball containing the manifolds M and N, 
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Figure 1: The time-dependent vector field v : fix [0,1] — >• R n+1 generates 
the morph F v (p,t), which is the solution of the initial value problem dq/dt — 
v(q,t), q(0) =p. 



C^°(fi;R" +1 ) the space of all smooth functions from fi to K" +1 with compact 
support, and V k := H^' 2 (fi; R n+1 ) the closure of C c °°(fi; R n+1 ) in the Sobolev 
space W k ' 2 (n-,R n+1 ) (see [§]). 

The space V k is a Hilbert space with the inner product 

n+l „ 

(/,fl>v* = X] E / D a f t D a g,dx, 

i=l a,H<fc ° 

where / = (/i, • • • , / n +i) : fi — * R n , a = (ai, . . . , a n+ i) is a multi-index 
with nonnegative integer components, |a| = a\ + . . . + ct rl +i, and fl"/; = 
d^fi/dx" 1 ...dx"^ 1 is the corresponding weak partial derivative of We 
choose k E N large enough so that the Sobolev space ,2 (fi) is embedded into 
C r (fi) and r > 1. By the Sobolev Embedding Theorem (see [2J[9]), it suffices 
to choose k > {n + l)/2 + r + 1. 

Consider time-dependent vector fields v : fi x [0, 1] — > K n+1 on fi that belong 
to the Hilbert space 

H k = L 2 (0,l;V k ) (3) 

(see Fig. IT]). By an abuse of notation, we will sometimes write v(x,t) — v(t)(x) 
for v € 7^ and (x,t) G fi x [0,1]. A time-dependent vector field v : fi X [0, 1] — > 
W n+1 belongs to the Hilbert space Ti k if its norm := (/* 

is finite. The inner product on Ti k is defined by 

(v,w) H k = (v(-,t),w(-,t)) V kdt. 
Jo 

Every vector field v E H k generates a morph F v : M x [0, 1] — > from 
M to F V (M, 1) via the evolution equation 

^ =«(«,*)• (4) 
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More precisely, let r] v (t; t , x) be the evolution operator of equation Q; that is, 
for every t € [0,1] and x € the function t i— > ri v (t;t ,x) solves equation Q 
and satisfies the initial condition rf (to; to,x) — x. The morph F v is defined by 
F v (p, t) = rj v (t; 0,p) for all (p, t) G M x [0, 1]. By the properties of the evolution 
operator rf , which have been studied in [5] and [2U], the morph F v (p,t) is of 
class AT' ac (M, F V (M, 1)) (see lemmas |bTT] and |B~2"t . The time-one map of the 



evolution operator rf is defined to be 4> v (x) := rj v (l; Q,x) for all x S 0, and we 
define ip v = <f> v \ M - 

Let „4|, be the admissible set of all time-dependent vector fields in 7i k that 
generate morphs between the manifolds M and N and are bounded by a uniform 
positive constant P. In symbols, 

A k P = {v e H k : t/j" E Diff r (M, N) and < P}. (5) 

We will prove that for P sufficiently large, the admissible set A k P is nonempty 
and Ap is weakly closed in H k (see lemma 2.6 1. 



Let T r S (M) denote the set of all continuous tensor fields on M contravariant 
of order r and covariant of order s (also called type (r,s)). For a tensor field 
ttv S T° s (N) and a diffeomorphism /i : M — > iV, /i*tat denotes the pull-back of 
rjv to M. 

For each t S [0, 1] and u e the manifold M 11 -* := F V (M, t) is called an 
intermediate state of the morph F v between manifolds M and TV generated by 
the time-dependent vector field v. We endow this intermediate state with the 
Riemannian metric g" t inherited from its embedding in K" +1 and let IIJ 1 denote 
the corresponding second fundamental form. 

Definition 1.3. Let B\ and E>2 be nonnegative real numbers, F v the morph, 
and <p v the time-one map generated by the time-dependent vector field v £ 
A P C TL k via the evolution equation Q. Recall that tp v :— 4> v \m- The bending 
distortion energy of v is 

E{v) = E(v;B 1 ,B 2 ) = B 1 [ \\(ip v )*g N - g M \\ 2 u M 

+ B 2 [ HW-TILv-IImHV/ 
Jm 



and the morphing distortion energy of v is 
£{v)=£{v\B x ,B 2 ) = B\ ( ( 

JO J A 



\F v (;tYg v t -g M \\ 2 uj M dt 



B 2 [ f \\F v (-,t)*ir t -II M fu; M dt, 
Jo Jm 



10 Jm 

where || • || is the fiber norm on the tensor bundle T° 2 (Af) generated by the 
fiber inner product g*M®9*M- (Note: We will use the same notation for the fiber 
norm on the tensor bundle T° S (M) generated by the inner product ®l = ig*M-) 

We will prove that the functionals E and £ have minimizers in Ap. 
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Theorem 1.4. (i) If P > and k E N are sufficiently large, then each of the 
functionals E : Ap — > R+ and £ : _4 P — > K + has a minimizer in the admissible 
set A P . 

The detailed conditions on the constants P and k are formulated in theo- 
rem 

We note that each difFeomorphism ip v : M — ► iV generated by a time- 
dependent vector field i> S .Ap is isotopic, as a map from M to R™ +1 , to the 
inclusion map i : M — ► via the isotopy S .M r ' ac (M, iV). To mini- 

mize the distortion energy of diffeomorphisms from other isotopy classes, we 
replace the map ip v in the definition of the functional E by the difFeomorphism 
ip v o <p : M — + N, where <fi is a fixed difFeomorphism on M. The existence of 
minimizcrs of the functional E with the above adjustment guarantees the exis- 
tence of minimizers of the functionals $ and A defined in displays ([lj and Q 
in a restricted admissible set of all C 2 diffeomorphisms between the manifolds 
M and N, which, considered as maps from M to R n+ , are isotopic to a given 
map 4>:M^ W n+1 . 

Theorem 1.5. If P > and k £ N are sufficiently large, then for every <f> G 
DifF(M) both functionals $ and A defined in displays ([T| and ^ respectively 
have minimizers in the admissible set 

Bp,4> : = i h e DifF 2 (M, N) : h= i) v o <f> for some v € A k P }. 

In section [3] we construct an example of a minimal distortion difFeomorphism 
and morph between the unit circle S 1 and the circle with radius R > 1, in 
R 2 . 

While the construction of a minimizer of the functional E does not cause 
significant difficulties, finding a minimizer of the functional £ is a much more 
intricate process. Even after we restrict our attention to the family of morphs 
whose intermediate states are concentric circles, finding a minimal distortion 
morph requires delicate analysis, which is done in section [3] 

To find a morph H(p,t) = i>{t)p with V € Q+ := {4> G C 2 (0, 1) H C[0, 1] : 
0(0) = 1, 0(1) = i?, and is increasing}, which has minimal distortion among 
the morphs F € A4 3 ' ac (M, N) whose intermediate states are circles with in- 
creasing radii, we solve the optimization problem 

minimize J(ip) = J^(ip 2 - l) 2 dt + - l) 2 dt 

for TpeQ+ ^ ( 6 ) 

subject to (iA dt< A, 

where A > 0. The inequality constraint in optimization problem ^ is derived 
from the requirement that the vector fields on the set £1 C M 2 generated by the 
morph H must be bounded by a uniform constant. 

More precisely, let be the open ball in K 2 of radius R + 2 and let p : K 2 — > 
K 2 be a bump function such that p = 1 on the open ball B(0, R + 1), p = on 
fi c , and < p < 1. 
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Given P > 0, define 

A(P) := Hpido \\ w ), (nm P 2 . 

Theorem 1.6. If the constant A — A(P) > log 2 R, then there exists a unique 
minimal distortion morph H(jp,t) — ip(t)p, where ip € Q+, between the unit 
circle S 1 and the circle of radius R > 1 in M 2 , among all the morphs F G 
M 3 ' ac (M,N) of the form 

F(p,t) = cf>(t)p, 
that generate the time- dependent vector field 

v{x,t) = ^§-p{x)x, (i,t)6fix[0,l] 

such that \\v\\-h.5 < P. 

Moreover, the radial function ip G Q+ of the distortion minimal morph H is 
the unique solution of the optimization problem (J6j) and solves the initial value 
problem 

/ = ^ + W 2 -i) 2 + W-i) 2 , m 
\^(o) = i, [) 

where the pair of positive constants A and fi is the unique solution of the system 
of equations 

^ (8) 



1 s 



and 



1 

V\Ji 



ty/fl-t 


- (s 2 - 


-1) 2 + 


(s-1) 2 




f (s 2 - 


-1) 2 H 


-(s-1) 2 


s 



ds = A. (9) 



2 Bending and Morphing via Time-Dependent 
Vector Fields in R n+1 



In this section we prove theorem |1. 4 1 

We will show that the admissible set Ap is nonempty if P is sufficiently 
large. 

Lemma 2.1. Let M and N be manifolds as in definition \l . 1\ Let F be a smooth 
morph between the manifolds M and N and assume that Q, C R™ +1 is an open 
ball in R™ +1 containing the image F(M x [0, 1]) of the morph F. There exists 
Pq > such that the admissible set A P is nonempty whenever P > Pq and 
k > S±i +2. 
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Proof. The morph F £ M(M, N) defines the R n+1 valued function 

r)F 

v(y,t) = —([F(.,t)]-\y),t) 

on the compact subset Q = {(F(x,t),t) : (x,i) £ M x [0, 1]} of R n+1 x R. 

We will extend the function v to a smooth vector field w £ Tt k such that 
^ W {M) = TV. 

First, notice that the smooth map G : M x [0, 1] -> x R defined by 

G(x,t) = (F(x,t),t) is a proper map (M x [0,1] is compact) and an injective 
immersion, hence an embedding (see [I]). Therefore, Q — G(M x [0,1]) is a 
submanifold (with boundary) of K n+1 x K (see [12] )■ 

Next, notice that the map G\ : Q — > Af x [0,1] defined by G\(y,t) = 
([F(-, i) is the inverse of G. Because G is an immersion, hence a local 

diffeomorphism, the map G\ is smooth. Therefore, the map v : Q — > IR n+1 is 
smooth because it is the composition of two smooth maps Gi : Q — » M x [0, 1] 
and If : M x [0,1] -> 

The smooth function w : Q — > ]R n+1 can be extended locally. That is, for 
every (y, t) £ Q there exists an open set J7 C K™ +1 x K such that (y, t) £ U and a 
smooth function : {/ — > M n+1 such that Vi\unQ — v \unQ- This local extension 
property follows from a more general fact about smooth functions defined on 
submanifolds: Let S be an s-dimensional smooth submanifold (perhaps with 
boundary) of M. m and let / : S — > R be a smooth function. Then for every 
x 6 S there exists an open set W C K m containing x and a smooth function 
fx : W — > K such that /|v7ns = /i I wns- It is easy to construct a local extension 
of the function / using submanifold charts on S and the definition of a smooth 
function whose domain is a submanifold with boundary. The details are left to 
the reader. 

Therefore, the function v : Q — > R n+1 satisfies the conditions of the smooth 
Tietze extension theorem (see pQ) and can be extended from the closed set 
Q C R n+1 xtbya smooth map v : R n+1 xR^ R n+1 . 

Finally, define w(x, t) — p(x)v(x,t), where p : R n+1 — > R is a smooth bump 
function such that p = 1 on Q and p = on d£l, and set Po := \\w\\^k. □ 

From now on, we assume that the open set in IR n+1 is chosen as in the latter 
lemma and the constant P > is large enough so that the set Ap is not empty; 
the number k of weak derivatives satisfies the inequality k > (n+l)/2 + r + l, 
where r > 1. 

For each v £ A k P , the time-one map ip v : M —> N transforms the interior 
(respectively, the boundary) of the manifold M to the interior (respectively, the 
boundary) of the manifold N. The existence and the convergence properties 
of the evolution operators generated by vector fields v £ TL k via the evolution 
equation Q have been studied in [HI HO]- For convenience of the reader, we state 
some of these properties ( which will be useful in our proofs) in Appendix [B] 

As mentioned in section 1.3 every time-dependent vector field v £ A% C TL 



generates a morph F v between the manifolds M and N of class A4 r ' ac (M, N). 
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Let us recall the distortion energy functionals E : Ti k — > R + and £ : H k 
l + (see definition 1.3 1. 



One of the main ingredients in the proof of theorem |1.4| is the weak conti- 
nuity of the functionals E and £ . We will prove the weak continuity of more 
general auxiliary functionals, where the tensor fields tm and in the following 
definition will later be replaced by the first and the second fundamental forms 
on the manifolds M and N respectively. 

Definition 2.2. Let M and N be manifolds as in definition |1.1| For given 
continuous tensor fields tm and tx of type (0, s) on M and N respectively, the 
functional J : A k P — > R+ is defined by 



J(v)= [ \WTrN-T M \\ 2 UM- 

JM 



Let v £ Ap be a time-dependent vector field that generates a morph F v € 
M. r ' ac (M, N) from the manifold M to TV. Recall that the intermediate state 
at the time t S [0,1] of the morph F v is denoted by M v,t . The Riemannian 
metric and the second fundamental form on M v,t , which are associated with 
the embedding of M v,t into R ra+1 , are denoted by g% and II" respectively. The 
functionals I\ : A k P — * M. + and I2 : A k P — » M+ are given by 

h(v)= I' [ \\F v (;tyg v t -g M \\ 2 u M dt (10) 



and 



f [ \\F V (;tyiI v t -II M fGJ M dt. (11) 
JO JM 



Definition 2.3. Let X be an n-dimensional vector space equipped with the 
inner product gx- Let T° S (X) — <S>f =1 X*. For every v E X, we denote the 
norm of v with respect to the inner product gx by \v\ gx — gxivjv) 1 / 2 and the 
unit sphere by S gx = {v € X : \v\ gx — 1}. Define the norm on T° S (X) by 

= m ax \b( Vl , . . . ,v s )\. 
Another norm on T° S (X) is defined by (see [15]) 

ll&ll = ®Ug*xQ>)- 

Let M(X) C T" 2 (X) denote the metric space of all inner products on X 
with the metric d(g, g') = \\g — g'\\ gx ■ 

Note that if {ei, . . . , e„} is an orthonormal basis of (X, gx), then 

n 

\\b\\= £ 6( eil ,...,e ls ) 2 . (12) 

ii ,.. . ,i s — 1 
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Lemma 2.4. The function r) : T° S (X) x M(X) -> K denned &y = WPWg 

is continuous on its domain. 

The proof of the latter lemma is sketched in Appendix [C] for completeness. 

For a C 1 Riemannian manifold (S,gs), let |j • || be the fiber norm on the 
bundle of all continuous tensor fields on S of type (0, s) generated by the fiber 
inner product <8>f = i<7s. 

Lemma 2.5. If b is a continuous tensor field of type (0, s) on a C 1 Riemannian 
n-manifold (S,gs), then the function 

z >-> \\b(z)\\g s (z) (13) 

is continuous on S. Moreover, the norms \\ ■ \\ and \\ ■ \\g SM are uniformly 
equivalent; in fact, 

-L\\ b ( z )\\<\\b(z)\\ gs(z) <\\b(z)\\ (14) 
for all b e T° S (S) and zeS. 



Proof. The continuity of the function defined in display ( 13 1 follows immediately 



from lemma 2.4 and, the inequalities in display (14 1 can be easily derived from 
the definitions of the norms || • || and || • || ss ( z ) and formula (12 1. □ 



Recall that a sequence {V}^ C Ti k converges weakly to v £ TL k (in symbols, 
v l — v) as I — > oo if (v — v, w) — > as I — > oo for every w E 7i . We call v the 
weak limit of -jV}^- A set Q C Ti k is sequentially weakly closed if it contains 
the weak limit of every weakly convergent sequence {v l }'^ 1 C Q. 

Lemma 2.6. (i) The admissible set Ap is sequentially weakly closed in Tt k . 
(ii) Let b be a continuous tensor field of type (0, s) on the manifold N; and, 
for every w € Ap, let ip w denote the restriction to M of the time-one map of 
the evolution equation dq/dt = w(q,t). If a sequence {v l } ( ^ 1 C Ap converges 
weakly to v e Ap in 7i k , then 

Km \\(r'yb-(r)*b\\(p Q ) = 

I — >oo 

for every p Q € M . 

Proof, (i) Let {v l }f^ 1 C A k P and suppose that v l converges weakly to some 
v G 1~C k as / — ► oo. We will show that v € A k P . 



By lemma B.3 rf (t;to,x) — * r/ v (t;to,x) as I — > oo (in the Euclidean norm) 
for all t, t £ [0, 1] and x € ft. Thus, the time-one maps generated by v l and their 
inverses converge pointwise: <fi v (x) — * <^"(a;) and (0" )^ 1 (a;) — * (<fi v )~ 1 (x) as 
/ — > oo for all x € il. Because the manifolds M and N are compact, <fi v (M) C iV 
and ((/) t ') _1 (7V) C M. In view of these inclusions, the C r diffeomorphism of 
f2 restricted to M is a diffeomorphism, that is, ip v — 4> v \m € Diff r (M, iV). 
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By passing to the limit as I — > oo in the inequality < (v — v,v) + 

P\\v\\ H k, it follows that |H| W fc < P. Therefore v € A P , as required. 

(ii) For simplicity, let us assume that s — 2. Let (£/,£) be a chart on M at 
Pq. It suffices to show that 

B l := W l yb(X,Y)(p Q ) - (^yb(X,Y)( Po ) -> (15) 

as / — > oo for all smooth vector fields X, Y on U. 
Using the notation 

q l = ^'(PO), 
Q l y = OfVoff'^fe), and 

Qy = D4> v Yo{4?)-\ y ) 



for all y € TV, the quantity -B in expression ( pL5[ ) is recast in the form 

B l = b(q l )(Z l ql ,Q l ql )-b(q)(Z q ,Q q ) 

= b(q l )(Z qll Q l ql )-b(q l )(Z qh Q l ql ) (16) 

+ b(q l )(Z qh Q l ql )-b(q l )(Z qh Q ql ) (17) 

+ b(q l )(Z qh Q ql )-b(q)(Z q ,Q q ). (18) 

Using definition |2.3| and noting that the Riemannian metric g^ is inherited 
from we estimate difference ( p~6[ ) as follows: 

b(q l )(Z ql ,Q l ql )-b(q l )(Z qh Q l ql ) < \\b{q l )\\ m{qi) \Z l ql - Z ql \ m+1 \Q l ql \ un+1 . 

By lemma B.3 \Z l ql — Z g z | Rn +i — ► as Z — * oo, and |Q^|rh+i is uniformly 
bounded in I € N. By lemma |2.5| there exists a constant C > such that 



l|k(<?0llc/iv(9 ! ) — C ^ or au> ^ e ^- Therefore, difference (|16| ) converges to zero as 
/ — > oo. Similarly, it can be shown that difference (17) converges to zero as 
/ — > oo. Difference (18) converges to zero as I — » oo because z i— > b(z)(Z z , Q z ) is 
a continuous function on [/. Hence, i?' — > as I — > oo. □ 

We say that a functional / : Ti k — ► R is weakly continuous on 7i fc if I(w') — > 
/(u) whenever the sequence -jV}^ C Ti k converges weakly to v in 7i k . 

Recall that the inequality k > (n+l)/2 + r+l guarantees the embedding 
of the Sobolev space W^ 2 (fl,R n+1 ) into C r (Cl,W n+1 ), where r > 1. 

Lemma 2.7. Assume that the constant P > is large enough so that the set 
A P is not empty. Let the Junctionals J,Ii,l2 be defined as in definition 

(i) If k > (n + l)/2 + 3, then the junctionals J : A% — * R+ and I\ : A P 
are weakly continuous. 

(ii) If k > (n + l)/2+4, then the functional 1 2 ■ Ap — * M+ is weakly continuous. 



2.2 
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Proof. Let {v 1 }^ C A P and suppose that v l con verg es weakly to some v £ H k 
as £ — > oo (in symbols — v £ TL k ). By lemma 2.6 w e „4p and J(v), h(v), 
and /2(f) are well-defined, 
(i) We will show that lim^^ J(v l ) = J{v). 

Let G := ® g* M . For tensor fields a, b £ T° 2 (M) and every p e M, we 
have the equality 

\\\a\\ 2 (p)-\\b\\ 2 (p)\ = \G(a + b,a-b)(p)\ < \\a + b\\(p)\\a - b\\(p). 

By applying the Cauchy-Schwarz inequality, we obtain the inequality 

\j(v i )-j(v)\ < f \\(r l yTN+(r)*TN-2TM\\\\(r')*TN-(r)*TN\\u>M 



M 



< (/ \W l Tr N + {rTr N -2r M f^ M ) 



1/2 



M 



X ( / \W 1 Ytn-WYt n \?um) 



1/2 



M 



By lemma [2~6 



lim ||(^rr J v-(^)Viv|| 2 (p)=0 



(19) 
(20) 



for all p £ M. 



Let K > be the constant in display (65 1 of lemma B.3 By lemma 2.5 



and because the manifold N is compact, there exists a constant C > such 
that 1 1 tjv (2;) 1 1 SJV (^) < C for all z £ N. Using the equivalence of norms ( 14 1 , we 
estimate 

\W 1 Ytn\\(p) < n^\W l Yr N {p)\\ gM{p) 

< n^\\T N (r l (p))\\ gN{ ^ (p)) \Dr'( P )\ S 

< n s/2 CK s . (21) 

Using inequalities ( 19 1 and pi} , limit (20 ), and the Dominated Convergence 
Theorem, we conclude that J(ir)—» J(v) as I — > 00. 

Let us show that the functional 1\ is weakly continuous. By an estimate 
analogous to (19 1, it suffices to prove the following statements. 

(I) lip £ M and t £ [0, 1], then 



lim (•,*)*&" -F°{-,tYgn J (p)=0' 

I — >oc 



(II) There exists Si > such that 

ii^Mr^ii 2 ^)^ 

for all p £ M, t £ [0, 1], and I £ N. 
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Because all the Ricmannian metrics are inherited from M n+1 , whose standard 
inner product is denoted by (•,•), we have 

F v * {;t)*g? (p)(X, Y) - F v (-,t)* g v t {p){X, Y) = (D X F V ' (p, t)X, D X F V ' (p, t)Y) 

-(D x F"(p,t)X,D x F v (p,t)Y), 

for all p £ M, X, Y £ T p M, and t £ [0, 1], where D x denotes the derivative with 
respect to the spatial variable. The right-hand side of this equation converges 
to zero as I — ► oo by lemma B.3 This completes the proof of statement (I). 
By the same lemma and inequality ( 14 1 , for every p G M we have 

\\F*\.,t)*gi\\{p) < n\\F*'(.,tygf(p)\\ gM(p) 

< n\D x F v ' (p,t)\ 2 

< nK 2 . 

This inequality implies statement (II). 

(ii) We will show the weak continuity of the functional I 2 - By an estimate 
analogous to (191, it suffices to show two facts: 

(III) If p G M and t G [0, 1], then 

lim ||F^(-,t)*Ilf -F v (;t)*IF t \\ 2 (l>) = 0. 



(IV) There exists S 2 > such that 

\\F^(;tyiif \\ 2 (p)<S 2 
for all p e M, t £ [0, 1] and I G N. 



We will first prove statement (IV). 

Consider a morph F w generated by a time-dependent vector field w £ Ap. 



By definition of morphs of class A^ r ' ac (M, N) in section 1.3 the orientation 
of each intermediate manifold M™' 4 , where t £ [0,1], is induced by the C 2 
diffeomorphism F(-,t) : M —> Al w,t . Let Af w,t (z) denote the unit normal to 
the manifold M w,t at the point z £ M w ' 1 . We assume that for every positively 
oriented basis of T Z M W '\ the set of vectors {X u . . . , X n , N w '\z)} is 

positively oriented in K™ +1 . 

For p £ M, let (U,£) be a chart at p, choose two smooth vector fields X and 
Y on U, and let 7 : [0, 1] — > U be a C 1 curve at p such that 7(0) = X p . It is 
evident that the inner product 

(^' t (^(7( S ),t)),D XJ F-( 7 ( S ),t)F 7(s) ) = (22) 

for every t, s G [0, 1]. Let us recall that for every f G [0,1] the function x i— > 
F w ' t (x,t) is defined for all a; £ ft and denote its second derivative at x £ ft 
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by D%F(x,t). By differentiating expression (22 1 with respect to s at s = 0, we 
obtain the equality 

F»(.,t)*ir?(p)(X p ,Y p ) = (V D ^ MX M w ^F w (p,t)),D x F w (p,t)Y p ) 

= -{M w '\F w { V ,t)),DlF w { P ,t)[X P ,Y p ]), (23) 

where V denotes the standard Riemannian connection on M. n+1 (see [TTj). 
For every p G M, let W p , Q p G T p M be unit length vectors such that 

\\Fv'(;tyiit( P )\\ gMip) = \F» l (;tyiif( P )(w p ,Q p )\. 



Using inequality (14 1 and equation (23 1, we have the estimates 

\\F^(;tyir t l \\( P ) < n\\F*'(.,tyir t \ P )\\ gMip) 

= n\F v '(;tyilf(p)(W p ,Q p )\ 

= n\ '\F vl (p, *)) , l£ J*' (p; [W r pj Qp] 

< nK. 

This complete s the proof of statement (IV). 



By lemma 



B.3 



if a G {0, 1,2}, then the derivative D"F V (p,t) converges to 
oo in the Euclidean norm for every p G M and t G [0, 1]. 



D«F v (p,t) as I 

Taking into account equation ( p3| , we sec that statement (III) follows from the 
convergence 

J^' t (F ul (p,t))^AT' t {F v {p,t)) (24) 
as I — > oo in ]R n+1 for every p G M and t G [0,1]. □ 

Theorem 2.8. Assume that the constant P > is large enough so that the set 
Ap is not empty. 

(i) If k > (n+l)/2 + 3, then there exists a minimizer of the bending distortion 
energy functional E in the admissible set A P . 

(ii) If k > (n + l)/2+4, then there exists a minimizer of the morphing distortion 
energy functional £ in the admissible set Ap. 

Proof. Let {v l }'^ 1 C A k P be a minimizing sequence of E, that is 



lim E(v l 



inf E(w) > 0. 



By lemma 2.6 the set Ap is sequentially weakly closed and bounded. Therefore, 
there exists a weakly convergent subsequence {v lk } c ^ 1 with the weak limit v G 
Ap. The functional E is weakly continuous by lemma 2.7 Therefore, E{v) = 
mf we ^k E(w) and v is a minimizer of E. 

The existence of minimizers for the functional £ is proved in the same fashion. 

□ 
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Remark 2.9. Theorem 



2.8 implies the existence of minimizcrs of the functional 



A defined in display ([2 1 in the admissible set 

Bp := {h e Diff 2 (Af, N):h= V for some v £ A k P }. 

The set Bp, among other maps, contains smooth diffeomorphisms / : M — > 
N C M™ +1 that are homotopic to the inclusion map i : M — > M n+1 and generate 
time-dependent vector fields in A k P . 

To minimize the distortion energy of diffeomorphisms from other isotopy 
classes, we consider the family of maps {ip v o <\> £ Diff r (M, N) : v £ A P }, 
where tj> is a fixed diffeomorphism of M. Similarly, given a smooth isotopy 
G : [0, 1] X M — > M, we consider the family of morphs {F£ £ M r ' ao (M, N) : v £ 
where Fg(p,t) = F v (G(p,t),t) for all (p,i) e M x [0, 1], as candidates for 
minimal distortion morphs. The most interesting example of this generalization 
is, perhaps, the case where G(p, t) — <fi(p) for some fixed diffeomorphism </> : 
M — > N, so that the admissible isotopies are from the class of morphs F v (4>(p),t) 
generated by time-dependent vector fields in A P , where p £ M and t S [0,1]. 

The latter idea leads to the definition of the functionals 

E 4> (y)=E 4> {v;B l ,B 2 ) = B l [ \\(ip v ° <P)*9n ~ 9m\\ 2 u M 

Jm 

+ b 2 f ||(v>w)*nj V -n J1 ,|| i W 



and 



£ G (v) = £ G {v;B 1 ,B 2 ) = Bt [ [ \\F G (;t)* g v t - g M f 

JO JM 

+ B 2 f f \\F£{;tyil v t -ll M fu M dt, 
Jo Jm 



where B 1 and B 2 are nonnegative real numbers (cf. definition 1.3 1, (j> <= Diff(M) 
and G : M X [0, 1] — > M is an isotopy. . 

Theorem |2.8| can be easily generalized to show that for P > and k 6 N 
sufficiently large, both functionals and £g have minimizers in Ap for every 
diffeomorphism <j) ■ M — ► M and isotopy G : M X [0, 1] — ► M. 

Theorem 2.10. Assume that the constant P > is large enough so that the set 
A P is not empty. Let (f> £ Diff (M) and let G : M x [0, 1] — > M &e an isotopy. 

(i) If k > (n + l)/2 + 3, iften iftere exists a minimizer of the bending distortion 
energy functional E^ in the admissible set Ap. 

(ii) If k > (n + l)/2+4, £/ien i/iere exists a minimizer of the morphing distortion 
energy functional £q i n the admissible set A P . 

(Hi) If k > (n + l)/2 + 3, iften 6o£/i functionals $ and A defined in displays (jl| 
and |2| respectively have minimizers in the admissible set 

Bp, 4, ■= V 1 e Diff 2 (M, N) : h — ijj v o 4> for some v £ A%}. 
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The latter theorem is an easy generalization of theorem 2.8 More precisely, 
let {b l }^ 1 be a sequence of tensor fields in T° S (M) such that lim^oo ||&'||(p) = 
and || 6 (p)|| < K for all p G M and I G N, where if is a positive constant and let 
4> G Diff(M). Then lim^oo ||0*6'||(p) = and there exists a constant K 1 > 
such that ||0*£r(p)|| < Ki for all p G M and Z G N. Using the above observation, 
lemma |2~7| is easily generalized to the case where ip v and F v are replaced with 
ip v o and respectively, and the proof of theorem 2.8 remains the same. 

the statement (iii), which is equivalent to theorem 1.5 



In theorem 2.10 



follows from the statement (i). 



3 A Minimal Distortion Morph 

We have proved the existence of minimizers of the functionals E and £, which 
produce minimal distortion diffeomorphisms and morphs between manifolds M 
and N. In this section, we consider the special case where M — S 1 is the unit 
circle in the plane and N = Sjj is the concentric circle of radius R > 1 and 
construct a minimal distortion diffcomorphism and morph between them. 

Our example of a minimal distortion morph in subsection |3.2| demonstrates 
the importance of the bound H^H^ic < P in the definition of the admissible set 
Ap. If this bound is not imposed, there is a minimizing sequence of morphs 
{i*n}5£Li such that the distortion energy 

*(F n ) := / / \\F n (;tyg?-g M \\ 2 UMdt 

JO JM 

+ 11 \\F n {;t)*lYl-ll M fLO M dt (25) 
JO JM 

tends to zero, where <?" and II™ are the first and the second fundamental forms 
of the intermediate manifold F n (M,t) induced by its embedding into E 2 . An 
example of such a sequence is F n (p,t) = <j) n (t)p for all t £ [0, 1] and p G M, 
where <j) n G C°°(0, 1) n C[0, 1] is a function whose values remain in the segment 
[1, R] and such that (j> n (t) = 1 for all t G [0, 1 — 1/n] and n (l) = R- From the 



representation of ( 25 1 in local coordinates (see ( 28 1 ) we derive 



9(F n ) = 2tt f [(cj>l - l) 2 + (</,„ - l) 2 ] dt (26) 
Jo 

< 27r[(i? 2 -l) 2 + (i?-l) 2 ]-; 

n 

hence, lirrin^oo ^(F n ) — 0. On the other hand, there is no morph H in the 
space M r > ac (M, N) with r > 1 such that = 0: otherwise, H(-,l) would 

be an isometry between M and N. The sequence {F n } ( ^L 1 converges pointwise 
to the discontinuous morph 



p, if < t < 1, 
Rp, if * = 1 
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whose distortion energy ^f(F) vanishes. 



Theorem 2.8 implies that every sequence of time-dependent vector fields 
{v n }'^' =1 C TL k such that each v n £ TL k generates the morph F n must be un- 
bounded in H k . 

In our example of a minimal distortion morph, we solve the optimization 
problem for the minimal distortion morph between S 1 and §> R in the class of 
morphs, whose intermediate states are circles of increasing radii, that are gener- 
ated by time-dependent vector fields whose norms are uniformly bounded by a 
positive constant P. Numerical solutions suggest that the second time-derivative 
d 2 F/dt 2 of the minimal morph F increases as P increases. In effect, the choice 
of the constant P in the definition of the admissible set A P sets a restriction on 
the magnitude of the curvature of the curves t <— ► F(p, t), where p £ M. 

We begin the construction of the minimal distortion morph with the example 
of a minimal distortion diffeomorphism between S 1 and E> R . This example is 
based on the theory of minimal deformation (as measured by the functional <&) 
bending of regular simple closed curves developed in [3]. 

3.1 A Minimal Distortion Diffeomorphism between two 
Circles 

We will construct a minimal distortion diffeomorphism between the circles M = 
S 1 and N = S R . 

For r > 1, we consider the functional A : Diff r (Af, N) — ► R+ defined in 
display Q2p . Also, using the radius R > 1 of 3 R , we define the radial map 
h R : R 2 -> R 2 by h R (p) = Rp. 

Lemma 3.1. The restriction of the radial map hp to S 1 minimizes the func- 
tional A : Diff 2 (M, N) -> R+ defined in display Q. 

Proof. Fix p G M and q e N, and let 7 : [0, L(M)) -> M and £, : [0, L(N)) -> N 
be the (positive orientation) arc length parametrizations of M and N respec- 
tively such that 7(0) = p and £(0) = q. The distortion energy functional A can 
be recast in the form 

f L{M) r L(M) 1 

A(u) = / (u 2 ~l) 2 dt+ / (-M 2 -l) 2 di 
Jo Jo R 

=: Ji(«) + J 2 (u), (27) 

where it = £~ 1 ohoj ; [0, L(M)) — > [0, L(N)) is a local coordinate representation 
of h £ Diff (M, N) with h(p) = q. By lemma 4.1 in [3], the functions ui(t) = 
L(N)/L(M)t and u 2 (t) = -i(A^)/L(M)t + L(A r ) minimize the functional Ji in 
the admissible set 

B= {u£ C 2 (0,L(M))r\C([0,L(M)]) : u is a bijection onto [0,i(iV)]}. 

The proof of the statement that u\ and 1*2 minimize the functional J2 in B 
follows along the same lines. 
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Therefore, the map hji\§i minimizes the functional ([2]) over the set of all 
maps h £ Diff 2 (M, N) such that, for our fixed p £ M, h(p) = Rp. 

If h £ Diff 2 (A/, N) is such that h{p) = q Rp, consider an isometry / : 
N — > N such that f(q) = Rp. Because f*gN = 9n and /* II jv = IIjv, we obtain 
A(/i) = A(/ o h) > A(/ir|§i), which proves the lemma. □ 

As before, let ip v denote the time-one map of the vector field v £ A k P re- 
stricted to M. Using lemma [XT] it is easy to construct a time-dependent vector 
field that minimizes the functional E(v) = A(ip v ) in the admissible set A k P . In 
fact, every vector field v° £ A k P that generates the time-one map <j> v such that 
its restriction to M is ip v — H p \m, minimizes the functional E. An example of 
such a vector field is 

v(x,t) = p(x)w(x,t) 
for all x in the open ball := 73(0, JJ + 2)cl 2 and t £ [0, 1], where 

w ur.t = ; —x 

y ' ' l + (R-l)t 

and p : R 2 — > M. is a bump function such that 0<p<l,p=lin the open ball 
B(0,R + 1) C M 2 , and p = on il c . The vector field u generates the morph 
F v (p, t) = (1+ (R— l)t)p, whose time-one map restricted to M is tp v (p) = Rp. 



3.2 A Minimal Distortion Morph 
between two Circles 

Let us assume, as before, that M = S 1 , N = S}j, and R > 1. 

In the previous subsection, we have constructed a minimizer of the functional 
E; the construction was quite straight-forward. The time-integral involved in 
the definition of the functional £ makes the construction of its minimizer a 
much more intricate process. We will restrict our attention to morphs that op- 
erate through images that are concentric circles, while leaving open the question 
whether a minimizer must be purely radial, as the problem of constructing a 
minimal morph within this family is difficult enough. Note that although our 
functional is formally defined in terms of time-dependent vector fields, it is the 
resulting morphs we will be working with directly. 

We will construct a minimal distortion morph between the circles M = S 1 
and iV = §)j in case R > 1. As before, let A / l 3,ac (M, TV) be the class of morphs 
between the manifolds M and N that are absolutely continuous in time and 
class C 3 in the spatial variable. Recall that for a morph F £ Al 3 ' ac (M, N) we 
define /* = F(-,t) £ Diff 3 (M, M*). We assume that the morph F is generated 
by a time-dependent vector field v £ A k P . 

Consider the functional ^ : M r ' ac (M, N) — > R + , where r > 1, defined by 

*(i ? )= / / ||(/ t )*g f - 5 o|| 2 ^A/^+ I I \\(f)*II t -II \\ 2 u; M dt, 
Jo Jm Jo Jm 
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where gt and lit are the first and the second fundamental forms on the inter- 
mediate state M* induced by its isometric embedding into R 2 . We notice that 



£{v; 1, 1) = $(F V ) for all v € A% (see definition |1.3[). 

Fix a point p E M. Let 7 be an arc-length parametrization of M that 
induces the positive orientation on M with 7(0) = p. Let be the arc length 
reparametrization of M 1 obtained from the parametrization /* o 7 such that 
£*(0) = /*°7(p) an d both £* and /*°7 induce the same orientation of M*. Such a 
parametrization can be obtained by solving the equation s(t, x) = y for x, where 
s(t, x) = Jq I/* o 7(r)| dr is the arc length function of the curve M*. Using the 
implicit solution x(t, y) of s(t, x) — y, we define £*(y) = /* 070 a:(t, y). Bec ause 
the morph F is generated by a time-dependent vector field i> € A P , lemma B.2 
implies that the function t 1— » Df t (p), where p 6 M, is absolutely continuous. 
It follows that the function t 1— ► £*(s) is continuous for every s € [0, L(M t )). 

The local representation of /* is given by u*(s) = (£*) _1 o/'o 7(s), where 
s e [0, 27r], and the energy of the morph F is 

2 " vd«S 2 



Jo 



nr^&y- 1 )'** <28> 



where : [0, L(M')] — > K is the curvature function of the intermediate state 
M t . 

Let us restrict our attention to the morphs whose intermediate states are 
circles of increasing radii such that each intermediate state M 1 of such a morph 
F is a circle of radius with ip S C 2 (0,1) n C[0, 1] a (strictly) increasing 
function. In symbols, 

ip S g+ :={<?(> S C 2 (0, 1) n C[0, 1] : 0(0) = 1, 0(1) = R, and is increasing}. 



The curvature function of M* is given by = l/ip(t). By lemma 3.1 the radial 
map between the circles M and M* minimizes the functional 



/*"/ \\(f t T9t-9o\\ 2 u M + IK/THt-IIoll^M. 

JM JM 



Therefore, 



V(F) > V{H) = 2ir [ - l) 2 dt + 2tt [ (tp - l) 2 dt, 
Jo Jo 

where the morph H e M oc,2 (M, N) is given by H(p, t) = ip(t)p. 

To determine the morph H(p, t) = ip(t)p of smallest distortion energy ^(H), 
we will minimize the functional ,/ : L 4 (0, 1) — * M+ defined by 



1 r i 

2 i\2 j. , / /„/, i\2 



J(ij>):= / {^-iydt+ / (V>-l) 2 dt (29) 
Jo Jo 
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over all admissible radius functions ip- To define the admissible set for the 
functional J, let us put this example into the context of time-dependent vector 
fields. 

Let f2 be the open ball of radius R + 2 in M 2 . Given a morph H (p, t) = ijj(t)p 
(where tp G Q + , t e [0, 1], and p € M), let us construct a time-dependent vector 
field v e H 5 = L 2 (0, l;W 5 ' 2 (fi;R 2 )) that generates H, where the number of 



weak derivatives k = 5 is chosen in view of condition (ii) of theorem 2.8 

Consider the class of morphs of the plane M 2 that have the form F(x, t) = 
ip(t)x, where tp £ Q+. Define a time-dependent vector field v : M 2 x [0, 1] — > M 2 
by 

dF 

v(F(x,t),t) — —(t,x) 

or, equivalcntly, 

-/ \ i>'(t) 

Clearly, the morph F satisfies the differential equation dq/dt — v(q,t). To 
obtain the required vector field v, multiply v by a bump function p : K 2 — > M 2 
such that p = 1 on the ball B(0,R+1), p = on fi c , and < p < 1. The vector 
field 

V{X,t) = -^-p(a;)a; 

belongs to the Hilbert space Tt k and generates the morph 

H(p,t) :=i>(t)p = F\ Mx[0}l] (p,t) 
for all (p,t) G M x [0,1]. 



In theorem 2.8 we require the admissible set Ap, for some fixed P > 0, to 
contain all vector fields u G TL k such that the norm of v is bounded by P and u 
generates a morph between the manifolds M and N. 

Therefore, in addition to the assumption that tp € Q + , we must assume that 
the time-dependent vector fields of the form v(x,t) = ^jj^p(x)x are bounded 
in Ti. k by a fixed constant P > 0. In symbols, the required bound is 

' * dt<p 2 . 



ll«ll^ = l|p-idn|| 2 <2(n;K2) ^ (^) 
After introducing the constant 



^• lda ll^ 5 ' 2 (n;R2) 



we obtain the constraint 



Girl := / (^-^j 2 dt — A< 0. (31) 
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Note that the functional J can be written in the form 



J(^)= / uty)dt, 







where the smooth function u(s) = (s 2 — l) 2 + (s — l) 2 is strictly increasing on 
(l,oo). 

To find a morph H(p, t) — ip(t)p with ip S Q+, which has minimal distortion 
among the morphs F £ .M 3 ' ac (M, AT) whose intermediate states are circles with 
increasing radii, we must solve the optimization problem 

minimize J(ip) 

for^eg+-{0eC 2 (O,l)nC[O,l] : 

0(0) = 1, 0(1) — R, and is increasing} ^ ' 

subject to G{4>) < 0. 



The solution of problem (32 1 is obtained using the following outline: We will 



consider the related optimization problem 
minimize J(ip), 

ip € Q 12 := {0 € W^ 2 (0, 1) : 0(0) = 1, 0(1) = R} (33) 
subject to G(0) < 0, 

where (because every function ip £ Q 1,2 is absolutely continuous on [0,1]) the 
boundary conditions in the definition of the set Q 1,2 are to be understood in the 
classical sense. We will determine the unique minimizer ip of the optimization 



problem (33 1 and show that -0 is an increasing C 2 function. Because Q+ C Q 1 ' 



the same function ip is the unique solution of optimization problem ( 32 1 



Lemma 3.2. There exists a unique solution ip £ Q 12 of the optimization prob- 



lem (331. Moreover, 1 < ip(t) < R for all t £ [0, 1]. 



Lemma |3.2| is proved using the direct method of the calculus of variations. 
First, we prove the existence of a minimizer of the functional J subject to the 
constraint G(ip) < in the admissible set W 1,4 / 3 (0, 1) with the appropriate 
boundary conditions, and then we show that the minimizer is, in fact, in class 
W^ 1,2 (0, 1). The inequalities 1 < ip and ip < R are proved by contradiction using 
the the cut-off functions h\{£) — max{l,^>(t)} and h,2{t) = mm{R,ip(t)}, which 
would yield smaller values of the functional J than the minimizer. The details 
are given in Appendix |D| 



Lemma 3.3. If the constant A in definition (31) satisfies the inequality A > 



(log R) (see also equation (30 \) andip £ Q is the solution of the optimization 



problem (33 I, then there exists a constant A > such that 



(i) ip is a critical point of the functional J + AG over the space of variations 
W 1,2 (0,1), and 

(ii) G{ip)={). 
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Moreover, the solution ip of the optimization problem (33 1 is in class C 2 (0, 1). 



Lemma 3.3 follows from the Generalized Kuhn- Tucker theorem (see theo- 
rem 1, Sec. 9.4 in [IB]) and a regularity result for weak solutions of Euler- 
Lagrange equations (see theorem 1.2.3, Sec. 1.2 in [13]). The proof of the lemma 
is sketched in Appendix [D] for completeness. 



Theorem 3.4. If the constant A in definition (31 \ satisfies the inequality 
A > (logi?) 2 , then there exists a unique function tp S C 2 (0,1) n Q 1 ' 2 satis- 
fying conditions (i) and (ii) of lemma 3.3 and the following properties, 
(in) The function ip is strictly increasing and solves the initial value problem 

(34) 



V/ = ^t/Vm + ^-i^ + W'-i) 2 , 

V(0) = 1, 



where the pair of positive constants A and fi is the unique solution of the system 
of equations 

71 (35) 



and 



1 

71 



ds 




(s 2 - l) 2 + 


( S -l) 2 




- (s 2 - l) 2 H 


-( S -l) 2 


s 



ds = A. 



(36) 



(iv) The function ip is the unique solution of the optimization problem (32 I. 

Proof. If ip e Z := C 2 (0, 1) n Q 1,2 is a critical point of the functional J\ := 
J + AG : W 1,2 (f), 1) — * M. + , then i/j satisfies the Euler-Lagrange equation for J\ 
which is equivalent to the Hamiltonian system 

f v' = f(<M, 

with the Hamiltonian H(tp,p) = pip' — L(ip, ip'), where 

m^') = (v/-i) 2 + (^-i) 2 + a(^) 



is the integrand of J\ and p := ■§^r(t/},ip > ) (see, for example, [7]). Moreover, 
the Hamiltonian H(ip,p) is constant along the solutions of the Euler-Lagrange 
equation for J\. Let us denote this constant by fi. 
It is easy to see that 



and the Hamiltonian is given by 

1 



H(iP,p) = -p 2 V - {V - \y - (V - iy 
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Note that the equation tp' — jj^(ip,p) yields ip' — ^pip 2 . By solving the 
Hamiltonian energy equation 



1 

4A 



for p and substituting, we obtain a first-order differential equation for 0: 

V/ = A=r/>y/» + W> 2 - I) 2 + (4> - I) 2 • 
v A 



(37) 



(38) 



The case with the negative square root is eliminated because the conditions 
ip(0) = 1 and '0(1) = R > 1 can be used to show that the derivative of ip is non 
negative on (0, 1). 



In view of equation (37 1, it is easy to see that fj, + (ip 2 — l) 2 + (ip — l) 2 > 
for all ip £ Z. Because 0(0) = 1, we have /i > 0. Also, it follows immediately 
from equation (38 1 that tp is an increasing function. 

Let us use the notation u(s) = (s 2 — l) 2 + (s — l) 2 and recall that u is a strictly 
increasing function on (l,oo). After integrating both sides of equation (38 1 over 
the interval < t < 1 and making the substitution s = ip(t), we obtain the 
relation 

" * (39) 
s^//i + m(s) VA 



Another relation of A and fi is obtained from condition (ii) in lemma 3.3 
(see equation (31 1 for the definition of G). The integrand in the definition 
of G contains the quantity (0') 2 , which we view as ip'i/j'. We substitute the 



right-hand side of equation ( 38 ) for one factor ip' of this square and leave the 



other factor tj/ in the resulting integrand. After making the change of variables 
s — ip(t), we obtain the equivalent relation 



1 

71 



^H + u(s) 



ds 



(40) 



We claim that there exists a unique solution (/i, A) of the equations (39 1 



and (40 1. To prove this, substitute for 1/vA from equation (39) into equa- 



tion (40 1 to obtain the equation 



a = m 



R y/fl + U{ S ) 



1 s 



ds. 



(41) 



Make the change of variables t = u(s) in both integrals in display (41 1 and 
then write /(/i) as a double integral to obtain the formula 



/(**) = 



«(«) MR) 



o Jo VJ^+~sH{t)H{s) 



ds dt, 



(42) 



where H(t) := it" 1 (t)u' (u' 1 (t)) > for all t £ [0,u(R)]. 
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By inspection of equation (41 1, it is easy to see that lim M ^ + /(/•*) = +°° 
and lim^-Kx, /(/i) = log 2 (i?). We will show that / is a decreasing function, 
which guarantees the existence of a unique solution of the equation /(/x) = A 
for all A > log 2 (i?). 

Using formula (|42|), we compute 



/V) = 2 



(fi + s^^ + t)^ H(t)H(s) 



ds dt. 



Let £>+ = {(s,t) € [0,w(i?)] 2 : s > t} and D_ = {(s,t) € [0,w(i?)] 2 : s < t}. 
After making a change of variables j(s, t) — (t, s), we see that 

s-t 1 

D+ (,/ + 3)3/2^ + f)l/2 ' 



-(is (it 



1 



_(/x + t) 3 / 2 (/x + s)V2 H(t)H(s) 



ds dt. 



Therefore, 
2/'(/i) 



t 



1 



( S -i) 2 1 



c?s (it 



D _ ( M + 5)3/2(^ + ^3/2 H(t)H( S ) 



ds dt < 0. 



This completes the proof that / is a decreasing function. 

There exists a unique solution /x of the equation /(/x) = A provided that 
A > log 2 (i?). The constant A is then easily found from equation (39 1. 

Having found the unique solution (/i, A) of the system ([39]) and (40 1, we solve 
the initial value problem (34). In fact, this initial value problem is equivalent 
to the integral equation 



</>(*) 



ds 



+ u(s) V A 



=t. 



(43) 



It follows that the unique solution tfj of the initial value problem (34 1 exists for 



all t £ [0, 1] and, because of condition (39 1, satisfies "0(1) = R- 



□ 



Figs. [2] and [3] depict graphs of the minimizer ip of the optimization prob- 
lem ^ with R = 2 and A = /(/x) in case /x = 0.001 for Fig. [2] and fi = 500 
for Fig. [3] Because / is a decreasing function of /x, Fig. [2] corresponds to a 
larger constant A. These plots illustrate that second derivative of the radius 
function ip corresponding to the minimal morph increases as the constant A in 
definition (31 1 increases. 
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m 




Figure 2: Graph of the radius function ip with R = 2, ^ = 0.001, A = 0.306067, 
and A = 1.56296. 




Figure 3: Graph of the radius function ip with R = 2, fi = 500, A = 1045.58, 
and A = 0.480456. 
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4 Minimal Deformation Bending 
of Two-Dimensional Spheres 

In this section we minimize the deformation energy functional <!> defined in 
display ([!]) under the assumptions M = S 2 and N = ft^(§ 2 ) =: RS 2 for some 
R > 0, where S 2 is the unit 2-dimensional sphere in IR 3 , and Hr is the radial 
map defined by hn(y) = Ry for all y £ IR 3 . As usual, the manifolds M and N 
are equipped with Riemannian metrics gu and gjq respectively induced by the 
Euclidean metric dx\ + dx\ + dx 2 of IR 3 . The manifolds M and N are Riemann 
surfaces (see p~7l [14] ) . We parametrize the spheres S 2 and R S 2 on the extended 
complex plane C = C U oo using stereographic projections. For (yi , j/2 , 2/3) £ S 2 , 
the stereographic projection is given by the expression 7r(t/x> Vii 2/3) — V i- l y 3 2 ■ 
We will show that maps of the form h = / oft H | M , where / is an isometry on AT, 
minimize the functional $ defined by equation ([l]) in the class of all holomorphic 
diffeomorphisms from M to N. 

We note that the holomorphic minimizers are critical points of the functional 
<f> on its natural domain Diff(M, N) (see appendix |A| and corollary A.4|. 

The parametrization <f> : C — > § 2 is given by 

2u 2v -l + u 2 + v 2 ^ T 



4>(u + iv)=[- 7, oi--, b-\ 1 ( 44 ) 

v ; V 1 + u 2 + v 2 ' 1 + u 2 + v 2 1 + u 2 + v 2 J y ' 

and the parametrization 4>r : C — > RS 2 of RS 2 is given by 4>r{u + iv) 

II'.')' :! + iv). 

In these coordinates, the Riemannian metrics gu and gjq are defined by 

4 

(1 + 

and 

4i? 2 

g N (z, z) = - dzdz. (46) 

Let ft, 6 Diff (M, AT) be a holomorphic map. The local representation (^_r) _1 o 
ho <j) : C C of h, which (by an abuse of notation) we shall denote by the same 
letter, is a holomorphic diffcomorphism of the extended complex plane onto 
itself. We conclude that h(z) has the form h(z) — M(z), where M(z) — [ s 
a Mobius transformation and a,b,c,d £ C are such that ad — be ^ 0. For such 
an h, it is easy to derive the formula 

u* , —\ AR 2 \bc-ad\ 2 ^_ 

h*g N (z,z) = - -^dzdz. (47) 

(\az + b\ 2 + \cz + d\ 2 ) 

Hence, the problem of minimization of the deformation energy functional $ 
defined in display ([T]) over all holomorphic diffeomorphisms from S 2 to RS 2 
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reduces to the problem of minimization of the function 



*(a, b, c,d) 



R 2 \bc-ad\ 2 1 \2 

2-7^ — ri9v? (! + k ) 



(48) 

where 2 = u + iv, over the group Aut(C) = PGL(2, C). Recall that the elements 
of the projective general linear group PGL(2,C) are the equivalence classes 
[a,b,c,d], where ad — be 7^ and {a! ,V ,d ,d!) € [a, b, c, d] if (a' ,b' ,d ,d') = 
X(a, b, c, d) for some A 6 C\{0}. 

Recall that the group of all isometries of the Riemann sphere is the projective 
unitary group PU(2,C); that is, every isometry / of (S 2 ,(?m) has the local 
representation (via stereographic projection) 



/(*) 



az 



cz + a 



where a,c € C are such that |a| 2 + |c| 2 = 1. 

The functional <!> is invariant with respect to left compositions with isome- 
tries; that is, $(/ o h) = $(h) for every isometry / € Diff(TV) and h 6 
Diff(M, N). Therefore, the reduced function ^ is well-defined on the quotient 
of PGL(2, C) by PU(2, C), which is the set of all equivalence classes 



a (3 
7 5 




a P 
7 8 



e PU(2,C) 



(49) 



We note that the equivalence class 





" 1 


" 









1 





consists of all the isometries of the unit sphere (S 2 ,(/m) 
Proof of theorem 11.21 



Statement (ii) of theorem 1.2 follows immediately from Hurwitz's automor- 
phisms theorem: The group of automorphisms of a compact Riemann surface 
of genus greater than one is finite (see |17j). 

Statement (i) of theorem 1.2 is equivalent to the following result. 



Theorem 4.1. The equivalence class of the isometries o/(S 2 ,<7m) is the unique 
minimizer of the function ^> defined on the homogeneous space 
PGL(2,C)/PU(2,C); that is, 



1 
1 



< * 



(50) 



for all 





a 


b ' 






c 


d 





E PGL(2,C)/PU(2,C). 
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Proof. The function W is well-defined on the homogeneous space PGL(2) / PU(2) . 
Thus, all values of if? arc obtained by choosing its domain to consist of one rep- 
resentative from each equivalence class. 
We claim that each equivalence class 



has a representative of the form 



a (3 
7 S 



1 

z r 



for some z E C and r E R + . 

To prove the claim, note that (without loss of generality) we may assume 
the determinant of the given representative is unity; that is, aS — [3j — 1. We 
wish to prove the existence of a, c £ C so that 



a (3 
7 5 



1 

z r 



(51) 



for some z E C and r E R + . In other words, it suffices to solve the system of 
linear equations 

aa — C7 = 1, 
a{3-c5 = 0. 



(52) 



In view of the equation aS — /3j = 1, it follows that a = S and c = (3. By 
substitution of a and c into equation (51 1, we find that z — (3a + and 
r = \(3\ 2 + \8\ 2 . This proves the claim. 

By the claim, it suffices to consider the value of ^ only at points of the form 
(1, 0, qe 1 ^, r), where q E R, r 6 R+, and ip E [0, 2ir). Thus, the theorem is an 
immediate consequence of the following proposition. 
The function * : R x [0, 2ir] x R + — > R given by 



*(g,V,r)=*(l,0,ge i *,r) 



(53) 



attains its global minimum on the set of points (0,ip, 1). 

To prove this result, let us first calculate the integral that represents the 
function 

After passing to polar coordinates (it = pcos 4> and v = psin </>), we represent 
"J in the form 



where £ = p 2 



oo />2tt 



^0 

pV- 



K 2 r 2 _JL 

(£ + f7COs(</>+'0)) 2 (1 + p 2 ) 5 



(1 + p 2 ) 2 pdct>dp, 



r 2 and r\ = 2pqr. Since the integrand is periodic with 
respect to (f> and we are integrating over one period, ty(q,tp,r) does not depend 
on ip; that is, 



V(q,>ip,r) 



OO />27T 



2^,2 



R 2 r 



1 



(£ + rj cos < 



(1 + p 2 ) 2 



(1 + p 2 ) 2 pdcj)dp. 



29 



The inner integral of the equivalent iterated integral is 



K(p) 



271 







R 2 r 2 



1 



(£ + ?7Cos(0)) 2 (1 + p 2 ) 2 

2tt 



(l + p 2 ) 2 pd0 



R 4 r 4 
- 2R 2 r 2 
f 2tt 



1 



o (£ + ?ycos(0)) 4 

2tt 



1 



1 



(l + p 2 ) 2 7 (£ + 77Cos(0)) 2 
' V(l + P 2 ) 2 - 



(1 + P 2 ) 4 



(54) 



Taking into account the inequalities £ > |r/| and > 0, the integrals in the 
previous expression are elementary; their values are given by 



1 



2rr£ 



{i + V cos{4>)Y 



and 



2tt 



1 



o (C + ??cos(0)) 4 



(£ 2 - r/ 2 ) 3 / 2 
<(2g 2 + 3r? 2 ) 

(£2 _ ^2)7/2 • 



By substitution into equation ( 54 1 , we find that 

r CO 

- / K(p)dp 



(55) 
./o 

= 7T - 2nR 2 + 7 ^ Y (l + q 2 + (r- l)r)(l + g 2 + r + r 2 ). 

The minimum of the function F(q,r) = ^(q^ip^r) on M x R + , is easily 
determined. Indeed, (0, 1) is the only critical point of F. Also, the Hessian of 
F is 

/ 47r_R 4 (l+3g 2 +r 2 ) 8TrR 4 q(l+q 2 ) \ 

D 2 F(qr) = 3r2 3r3 

W ' j I _8rfW) 2^fl 4 (3+6 ? 2 +3g 4 +r 4 ) J 

\ 3 r 3 3r 4 / 

We note that d 2 F(q,r) jdq 2 and the determinant of the Hessian 



det (D 2 F(q,r)) 



9r 6 



(3 + q 6 + 3r 2 + r 4 + r 6 + g 4 (5 + 3r 2 ) + g 2 (7 + 6r 2 + 3r 4 )) 



are both positive by inspection. By Sylvester's criterion, the Hessian is positive 
definite over the entire domain of F; therefore, F is convex. If follows that (0, 1) 
is the unique global minimizer of F. The minimum of F is 

F(0, 1) = ir(R 2 - l) 2 . 

Clearly, points of the form (0, tp, 1) are the global minima of on its domain 
R x [0,2tt] x E+. □ 
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5 Discussion 



For diffeomorphic hypersurfaces M and N in a Euclidean space, we have denned 
functionals that measure how well a diffeomorphism ip : M — > N preserves the 
geometry of M and proved that minimizers of these functionals exist. Since 
our functionals involve comparisons of the first and the second fundamental 
forms on M with the pull-backs of the corresponding fundamental forms from 
N to M by ip, they measure the extent to which ip changes the size of vectors 
carried from M to N and the extent to which it preserves the amount of bending 
of the unit normal vector field on M. In addition, to maintain flexibility for 
applications (where some particular combination of the measurements given by 
the first and the second fundamental forms is desired) we allow the measurement 
of distortion length and distortion bending to be weighted. 

Even in the case where M and N are one-dimensional, there are examples of 
such manifolds where the distortion energy functional for the first fundamental 
form has no minimizer in Diff (M,N) (see [3]). This fact and the complexity of 
the corresponding Euler-Lagrange equation (derived in appendix [A]) evidence 
the difficulty of the general problem of minimization of distortion energy func- 
tionals over Diff (M, N). 

One of the main ideas used successfully here is to restrict , via an appropri- 
ate linearization, the space of diffcomorphisms on which the distortion energy 
functionals are minimized. More precisely, instead of working directly with 
maps in Diff(M, N), we consider time-dependent vector fields on the ambient 
Euclidean space that generate the desired diffcomorphisms. Each such vector 
field determines an evolution family that can be applied to M, for values of the 
evolution parameter in the interval [0,1], to produce a morph carrying M to 
some diffeomorphic end-hypersurface. We restrict our attention to the class of 
vector fields for which this end-hypersurface is N. The advantage is clear: The 
elements of Diff (M, N) are replaced by elements of a vector space. Using this 
approach, which already appears in the literature on image deformation (see 
[Um[10j[19]), we are able to prove the existence of minimizers for the distortion 
morphing functional £ and the distortion bending functional E. The functional 
£ is an appropriate choice to measure distortion in case we wish to include the 
deviation of the geometry of the intermediate surfaces from the original surface 
M. On the other hand, the distortion functional E is appropriate if we wish to 
ignore the distortions represented by the intermediate surfaces and only consider 
the distortion caused by mapping M to N. 

We have also gone beyond the existence of minimizers. Indeed, we have 
determined minimizers of our distortion energy functionals for some classes of 
one-dimensional manifolds and for Ricmann surfaces. While these special cases 
have independent interest, they also serve to test conjectures. 

The existence of minimizers provides a hunting license for seeking approxi- 
mations of these minimal distortion diffeomorphisms and morphs via numerical 
methods. The known minimizers for special cases can be used to test numerical 
algorithms. 

Our admissible set ^4p of time-dependent vector fields is the closed ball of 



31 



radius P > in the Hilbert space TL k of vector fields intersected with the set 
B of all vector fields that carry M to N. The boundedness of A k P in H k is an 
essential ingredient in our proof of the existence of distortion energy minimizers. 
Thus, there is a natural (open) question: Is there a minimizer in the interior of 
A P (in the relative topology of the B) or are all minimizers on the boundary of 
this set? 

We mention some evidence that suggests there are energy minimizers of the 
functional E in the interior of A k P for P > sufficiently large. We expect the 
same result for £ , but this case seems to be much more subtle. 

Our result, discussed in section 3.1 that (1) the radial map is a minimal 
distortion map in Diff (S 1 , S^) in case R > 1 and (2) there is a time-dependent 
vector field whose time-one map is the radial map, shows that (at least for this 
special case) interior minimizers exist. Indeed, for P larger than the Hilbert- 
norm of our time-dependent vector field that produces the radial map, this 
vector field is in the interior of the admissible set A k P . 

To prove the existence of an interior minimizer in general, it suffices to show 
that there is a number Pq > such that for every time-dependent vector field 
v £ H k that carries M to N such that its Hilbert-norm IH^fe > Pq, there 
exists a vector field v' whose Hilbert-norm is less than Pq and whose energy (as 
measured by E or £) does not exceed the corresponding energy of v. 

Suppose that for some large P we have a minimizer of the functional E, a 
time-dependent vector field v whose Hilbert-norm is P; that is, v lies on the 
boundary of the admissible set. The time-one map of v determines the value 
of E. We seek a new vector field v' with smaller Hilbert-norm and the same 
time-one map as v; it would be our desired interior minimizer. A large Hilbert- 
norm for v suggests that its integral curves have large lengths or largeness in 
some measure of bending that takes into account the space derivatives of v 
which are used to compute the Sobolev norm. By shortening and straightening 
the integral curves of v in some subset of the space f2 x [0, 1], where our time- 
dependent vector fields are defined, we can construct a new time-dependent 
vector field v' whose Hilbert-norm is strictly smaller than the norm of v and 
whose time-one map is the same as the time-one map of v. 

Our choice of admissible time-dependent vector fields, which are defined 
on x [0, 1], for a ball in R n+1 containing the manifolds M and N, can be 
replaced by other sets of functions chosen to not contain extraneous information. 
For example, it might be desirable to consider only the curves leading from M 
to N and not all the integral curves generated by our time-dependent vector 
fields defined on f2 x [0,1]. One alternative admissible set is a subset of a 
class of functions we call development vector fields. They are defined to be 
functions v : M x [0, 1] — ► T1R™ +1 (in an appropriate function space) that 
generate morphs F v : M x [0, 1] — > ]R n+1 via the following construction. Let 

: T p W l+1 — > T q M. n+1 be the parallel transport from the tangent space at p to 
the tangent space at q defined by the Euclidean metric (or perhaps some other 
metric) on K n+1 . For each p £ M C the solution of the initial value 

problem 

c P = Tp p v(p,t), c p (0)=p 
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is a curve in M. n+1 starting on M at p. Thus, we have defined a morph F v : 
M x [0,1] -> R n+1 given by F v (p,t) := c p {t). The new admissible set is all 
development vector fields that are in some closed ball of radius P > of the 
corresponding Hilbert space such that F V (M, 1) = N. 

We also mention the possibility of defining a new norm on our time-dependent 
or development vector fields by reversing the order of integration: 



where 1Z is f2 for time-dependent vector fields and M for development vector 
fields. This norm may differ from the norm ||w||^fc used throughout this paper as 
functions in the Hilbert space Tt k are not necessarily L 1 in the joint (a;, t) variable 
on f2 x [0, 1]. Possible advantages of this approach are a better understanding 
of the relation between this norm and the shapes of the (integral) curves that 
are used to define morphs and, by first integrating over the temporal parameter, 
the Sobolev-norm can be viewed as a norm for functions defined on the compact 
manifold M. Consideration of this norm led us to the concept of a development 
vector field. Alternative approaches to the problem of minimal morphing using 
these ideas are a subject for future research. 

A The Euler-Lagrange Equation for the 
Deformation Energy Functional 

We will determine the Euler-Lagrange equation for the deformation energy func- 
tional $ : Diff (M, N) -> R + defined in display @. Let c : (-£, e) -> Diff (M, N) 
be a C 1 curve at ft. € Diff(M, N), which we call a variation of h. The equiva- 
lence class [c]h E Th Diff (Af, N) can be identified with the smooth vector field 
Y € X(N) defined by Y(q) = f t c{t) o h~ l (q) for all q e N. We call the vector 
field Y E X(N) a variational vector field of h € Diff (A/, AT). We will compute 
the first variation D$>(h)Y for all directions Y E X(N). 

Consider a smooth vector field X E X(Af) with flow 4> t and a diffeomor- 
phism h E Diff (M, N), and suppose that the variation c(t) = ho <p t induces the 
variational vector field Y = h*X E X(N). The diffeomorphism h is a critical 
point of the functional $ if and only if 



~$(hoct> t )\ t =o=D<S>(h)KX = 2 / G{h*g N -g M ,Lxh*g N )LJM = (56) 



for all X E X(M), where G = g* M ® g* M . 

Let V and V be Riemannian connections on M compatible with the Rieman- 
nian metrics a = gM and [3 = h* g^ respectively, and denote the corresponding 
Christoffel symbols of V and V by T l - k and t l - k . 

Let Y be a smooth vector field on M expressed in components by Y = 
Y k jf~k. The components of the Lie derivative of the Riemannian metric (3 in 
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the direction of the vector field Y are 

[Lyflkm = V fc K m + V m r fe , (57) 

where Y m — [0\ m jY 3 are the lowered coordinates of Y via the Riemannian 
metric (3 (see PQ). 

Recall that T r s (M) is the set of all continuous tensor fields on M contravari- 
ant of order r and covariant of order s, or type (r, s). 

Definition A.l. (i) Define the tensor field B = (f3 — a)** G T 2 (M). In 
other words, i3 equals the strain tensor field h* g^ — gu with both indices raised 
via the Riemannian metric a = gu- Its components are given by B km — 
OSy - n, ; in'''n-"<\ 

(ii) The bilinear form A : X(M) x £(M) -> X(M) is defined by 

A(x,y) = v x r-v x r (58) 

for all X,y G 3C(M). 

Remark A. 2. The bilinear form A can be viewed as a tensor field of type (1, 2) 
on M with components 

A m kp = VT p -TT p (59) 

(see [T5], proposition 7.10). 

Recall that the divergence of a tensor field r G T' S (M) is defined to be (see 

divr = C; +1 (Vr), (60) 

where C\ denotes the contraction in lower i and upper j index. The divergence 
of r, divr, is a tensor of type (r — 1, s). 

For two tensor fields G T r 2 {M) and r G T 2 S (M), 9 : t denotes the type 
(r, s) tensor field obtained by the contraction of the two covariant degrees of 
with the two contravariant degrees of r. 

Proposition A. 3. Let M be a compact, connected, and oriented smooth Rie- 
mannian n-manifold without boundary isometrically embedded into R n+1 . For 
the functional $>(h) = J M \\h*gjsr — gu || 2 % with domain Diff [M, N), we have 
that 



D$(h)(h*Y) = -4 / g M (divB + A:B,Y)cj M (61) 

JM 

for all vector fields Y G X(M), where the tensors A and B are as in defini- 
tion A.l. Moreover, h G Diff (M, TV) is a critical point of the functional $ if 



and only if 

divB + A : B = 0. 
The latter equation can be rewritten in components as follows: 

d k B km + T p kp B km + tt p B kp = (62) 

for all m = 1,2, ... ,n. 
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Proof. For given Y G X(M), consider the vector field 



X = (/3y - aij )a ik a jm Y m A, (63) 

where a — gu and /3 = h*g^. Although we describe X pointwise using local 
coordinates, X is a well defined smooth vector field on M because it is obtained 
by various contractions of the tensor fields a, (3, and Y. It can be verified that 
the divergence of the vector field X with respect to the Ricmannian metric gu 
is 

1 

2* 

Because M is without boundary, J M div gM Xojm = 0. Using this and the equal- 
ity 



dW gM X=-G(p-a,L Y 0) + g M (divB + A:B,Y). 



D$(h) (h*Y) = 2 / G{p-a,L Y /3)ujM, 

JM 



we conclude that 



D$(h)(h*Y) = -4 / g M {divB + A: B,Y)oj m 

JM 

as required. □ 



Corollary A. 4. Let M be a manifold as in proposition A. 3 and h £ Diff(A/, N). 
If h*9N = R 2 gM for some R G M, then h is a critical point of the functional 

In particular, let hjj : R™ +1 — > M n+1 &e the radial map defined by hn(p) = Rp 
for all p G where R > 0. Assume that N = hii(M) is a reseated version 

of the manifold M, and the Riemannian metrics g^ and gpj on the manifolds 
M and N are inherited from M. n+1 . Then every composition h = f o hji\M of 
the radial map }ir\m with an isometry f G Diff(TV) is a critical point of the 
functional <&. 



The proof of this corollary is an easy verification of equation ( 62 1 



B Existence and Convergence 

Results for Evolution Operators 

In this section we state results on existence and convergence of certain evolution 
operators. 

We denote the Euclidean norm of an element A G M m , where m G N, by 
\A\ and the Hilbert space L 2 (0,l;V k ) by H k , where the Sobolev space V k — 
Wa' 2 (n-W l+1 ) is embedded into C r (fi;E n+1 ) and r > 2. Recall that Sobolev's 
theorem guarantees the latter embedding if k > (n+l)/2 + r+l. The following 
lemma is proved in [5] . 
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Lemma B.l (Dupuis, Grenander, Miller). For every time- dependent vector 
field v G 7i and to G [0, 1], there exists a function tj> : [0, 1] X fi — ► W l+1 such 
that t i— ► </>(£, a;) is ifte unique absolutely continuous solution of the initial value 
problem 

g(*o) = x lD4j 

/or all t G [0, 1]. Moreover, the function x i— > cf>(t,x) is a homeomorphism of CI. 

For every d € 7i fe and x € 0, let F v (x,t) be the solution of the evolution 
equation dq/dt = v(q,t) with the initial condition F v (x,0) = x. For a function 
/ G C r (fl), denote 

\\f\\r,oo= ™p\D a f(x)\, 
cz,\ot\ <r 

where a = (ai, . . . , a n +i) is a- multi-index with nonnegative integer components, 
|a| = ai + . . . + a n+1 , and D a f = feQ1 J, n+1 . 

More general versions of the following two lemmas are proved in [20 , Ap- 
pendix C]. 

Lemma B.2 (Trouve, Younes). If v G H k and F v : Cl x [0, 1] -> R" is de/zned 
as above, then the function x i— > F v (x,t) is in class C r (Cl) and, for all q < r, 

~DlF*{x,t)=Dl(v{F v {x,t),t)), 

where denotes the derivative with respect to x of order q. Moreover, there 
exist positive constants C and C such that 



sup \\F v (;t)\\ rt00 <Ce c 'NI** 
te[o,i] 



for all v G 



Recall that we say v ; — ^ w weakly in Ti k as Z — > oo if («' — u, w) — > as 
Z — > oo for all w 6 7i fe . 

Lemma B.3 (Trouve, Younes). If the sequence {v l }f^ 1 C TL k converges weakly 
to v £ H k as I — > oo, iften 

sup ||F t,! (-,t)-F t '(-,t)|| r _ 1 , 00 ^0 
te[o,i] 

as Z — > oo. Moreover, there exists a constant K > smc/i £/iai 

sup ||F u '(-,t)|| r ,oo <if (65) 
te[o,i] 

/or a// Z G N. 



3G 



C Proof of Lemma 



2.4 



Proof of lemma 12.41 

Proof. The continuity of rj in the [3 variable is an immediate consequence of the 
equivalence of all the norms on the finite-dimensional space T° S (X). 
Let us show that rj is continuous in the g variable. 

Fix (3 G T° S (X) and g e M(X). Let K(g) be a positive constant such that 
\ v \g x — f° r an v € Sg, and let C(g) be a positive constant such that 

\v\g x < C(g)\v\ g for all v £ X. 

Choose e > 0. We leave it to the reader to show that there exists S > such 
that for all v G S g and g' e M(X) satisfying ||p — </|| Sx < <5 we have 

(66) 



IV S ||/3|| 9X C( 3 )(2X(. 9 )) 



s-1 ' 



Let the vectors Vi,...,v a £ S g be such that \(3(vi, . . . ,v s )\ = \\f3\\ g , and 
define m := Vi/\vi\ g i € 5 S ', where — <7'|| gx < <5- Then 

= !/?(«! -U X +U 1 ,...,V B -U s +U S )\ 



< . . .,u.)| + ||/5|| Sx C(.g)(2X( ff )) s - 1 £ ^ - Wi | 

< lulls' +£• 



»=i 



Therefore, ||/3|| ff ' > \\{3\\ g — s. The inequality \\(3\\ g > \\(3\\g> — £ can be shown 
in the same fashion. Hence, the function rj is continuous in the g variable. □ 



D Proofs of lemmas 3.2 and 3.3 



Proof of lemma 13.21 

Proof. The proof consists of two main steps: (1) Using the direct method of 
the calculus of variations, we will prove the existence of a minimizer for the 
auxiliary optimization problem 



minimize J(V0i 

tjj E Q M/3 := {(f>e W^i(0, 1) : ip(0) = 1,^(1) = R} 
subject to G{ij)) < 0. 



(67) 



(2) We will show that the minimizer for problem (67) is in W ' (0, 1) 



If follows that this minimizer is a minimizer of the optimization problem ( 33 ) . 

Let {ipn}^Li C Q 1 ' 4 / 3 be a minimizing sequence for the optimization prob- 
lem (67 1. In particular, G(tp n ) < for every positive integer n. In symbols, 



J(V>») - inf J(^). 

V-eQ 1 ' 4 / 3 ,C3(i/')<o 
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We claim that the minimizing sequence is bounded in W 1 ' 4 / 3 ^, 1). To prove 
this fact, we use the triangle inequality for the L 2 (0, 1) norm to make the esti- 
mate 

{[\uty /2 = {A^-i+i)w /2 

Jo Jo 

< J(V^„) 1/2 + i 

< Vm + i, 

where M > is a uniform bound for the convergent sequence {J(ip n )}^Li- By 
Holder's inequality with the conjugate constants 3 and 3/2, 



< 



< 



((Vm + i)a) 



2/3 



as required. 

Because the Banach space H Al ' 4//3 (0, 1) is reflexive, ip n — ip weakly in 
W 1A ^(0, 1) for some ^ £ H /1,4 / 3 (0, 1), up to a subsequence. We have ■0 £ Q 1 ' 4 / 3 
because the subspace Wq 3 (Q, 1) is weakly closed in W /1,4 / 3 (0, 1). 

The integrands {ip 2 — l) 2 + — l) 2 and {ip') 2 /ip 2 of J and G respectively 
are both convex functions of ip' '. Therefore, the functionals J and G are weakly 
lower semicontinuous in iy 1,4 / 3 (0,l) (see theorem 1, Sec. 8.2 in [3]). But then 



G{ip) < liminfn^oo G{ip n ) < and ip £ Q 1,4 / 3 solves optimization problem (67 1 
To prove that ip > 1, let us assume, on the contrary, that there exists (in the 
usual topology of [0, 1]) an open set W of positive measure such that ip{t) < 1 
for all t £ W. Define the cut-off function hi £ Q M/3 by h x {t) = max{l>(t)}. 
It is easy to check that G{h\) < and that J(hi) < J{ip), which contradicts 
the minimizing property of ip. The inequality ip(t) < R for all t £ [0, 1] can be 
verified in a similar fashion, using the cut-off function Ji2{t) = min{_R, ip{t)}. 
Using the inequality ip < R, we have the estimate 

~Wfdt = J^ 2 (j) 2 dt<R 2 A. 

Therefore, ip belongs to the space M^ 1,2 (0,1). 

Finally, the uniqueness of ip follows from the fact that the equality J{ip\) — 
J {^2)1 where tpi,ip2 G Q 1 ' 2 are such that 1 < ip 12 < R, implies uo ip\{t) — 
uoip 2 (t) for all t £ [0, 1], where the function u(s) — (s 2 — l) 2 + (s — l) 2 is strictly 
increasing on (l,+oo). 

□ 

Proof of lemma 13.31 
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Proof. Statements (i) and (ii) follow from the generalized Kuhn- Tucker theorem 
(see theorem 1, Sec. 9.4 in [IS]). 

We will verify that the minimizer ip is a regular point of the inequality 
G(ip) < 0. We leave it to the reader to verify that the functionals J and G are 
Gateaux differentiable at ip > 1 . 

It suffices to show that there exists h £ W ' (0, 1) such that 



o v> 3 



6G(i/>,h)= I -^(h'ip-ip'h)dt<0, 



where 5G(ip, h) is the Gateaux derivative of G in the direction h. 
Assume, on the contrary, that 



o ip 



3 (h'tp - ip'h) dt = (68) 



for all h e Wq 1,2 (0, 1). Then ip satisfies the Euler-Lagrange equation for the func- 
tional G whose associated Lagrangian (ip' /ip) 2 has a positive second derivative 
with respect to ip'. By a regularity result for weak solutions of Euler-Lagrange 
equations (see theorem 1.2.3, Sec. 1.2 in [13]), ip is of class C 2 (0,1). Therefore, 



we can integrate by parts in equation ( 68 1 to obtain the differential equation 

r = (^') 2 /V'- 

The function £ i — ^ J?* is the unique solution of the latter differential equation 
satisfying the boundary conditions ip(0) = 1 and — R- Therefore, the 



solution ip of the optimization problem (33 1 must be ip(t) — Rr. But, there is a 



function hp G Q ' such that G(hp) < and J(hp) < J(ip), in contradiction to 
the minimizing property of ip. In fact, a family of such functions is given by 

f 1, Xte[O,0\; 
M*):=S ^(t-P) + h if* S 09,2/?]; 
I R*, if t e (2/3,1] 

for f3 > sufficiently small. 

The C 2 regularity of the solution ip of the optimization problem ( 33 1 follows 
from (i) and the special form of the Lagrangian 



L(q,p) = (q 2 -l) 2 + ( q -l) 2 + \(P) 



2 



associated with the functional J+XG for A > 0: it has positive second derivative 
with respect to p on a neighborhood U of the set {(ip(t), ip'(t)) : t € [0, 1]} (see 
theorem 1.2.3, Sec. 1.2 in [H]). □ 
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